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Many experiments in the near future will test dark energy through its effects on the linear growth 
of matter perturbations. In this paper we discuss the constraints that future large-scale redshift 
surveys can put on three different parameterizations of the linear growth factor and how these 
constraints will help ruling out different classes of dark energy and modified gravity models. We 
show that a scale-independent bias can be estimated to a few percent per redshift slice by combining 
redshift distortions with power spectrum amplitude, without the need of an external estimation. We 
find that the growth rate can be constrained to within 2-4% for each /S.z = 0.2 redshift slice, while 
the equation of state w and the index 7 can be simultaneously estimated both to within 0.02. 
We also find that a constant dimensionless coupling between dark energy and dark matter can be 
constrained to be smaller than 0.14. 



I. INTRODUCTION 



The linear growth rate of matter perturbations is one of the most interesting observable quantities 
since it allows to explore the dynamical features related to the build-up of cosmic structures beyond 
the background expansion. For example it can be used to discriminate between cosmological models 
based on Einstein's gravity and alternative models like f{R) modifications of gravity (see e.g. [1]) or 
multi-dimensional scenarios like in the Dvali-Gabadaze-Porrati (DGP) [2] theory (e.g. [3] and references 
therein) . In addition, the growth rate is sensitive to dark energy clustering or to dark energy-dark matter 
interaction. For instance, in models with scalar-tensor couplings or in /(i?) theories the growth rate 
at early epochs can be larger than in ACDM models and can acquire a scale dependence [4-6] (see for 
instance [7] for a review on dark energy). 

Simultaneous information on geometry and growth rate can be obtained by measuring the galaxy power 
spectrum or the 2-point correlation function and their anisotropics observed in redshift space. These 
redshift distortions arise from peculiar velocities that contribute, together with the recession velocities, 
to the observed redshift. The net effect is to induce a radial anisotropy in galaxy clustering that can be 
measured from standard two-point statistics like the power spectrum or the correlation function [8]. The 
amplitude of the anisotropy is determined by the typical amplitude of peculiar velocities which, in linear 
theory, is set by the growth rate of perturbations^: 

s^P^. (1) 

where G{z) = 5{z)/5{0) is the growth function, 6{z) the matter density contrast and the scale factor a is 
related to the redshift z through a = (l-l-z)"-*^. Since however we only observe the clustering of galaxies 
and not that of the matter, the quantity that is accessible to observations is actually 

(2) 

where the bias b is the ratio of density fluctuations in galaxies and matter. The bias is in general a function 
of redshift and scale, but in the following we will consider it as a simple scale-independent function. 

Once the power spectrum is computed in fc-space, the analysis proposed in [9] can be exploited to 
constrain not only geometry but also the growth rate (as pointed out in [10]; see also ]3, 11]), provided 



^ In order to avoid confusion with the f{R) models, we use the letter s, for slope, rather than the more popular /, to 
indicate the growth rate. 
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that the power spectrum is not marginahzed over its amplitude. In configuration space, the first analysis 
of the two-point correlation function explicitly aimed at discriminating models of modified gravity from 
the standard ACDM scenario has been performed by [12]. Currently, there are several experimental 
estimates of the growth factor derived from the analysis of the redshift space distortions [12-19], from the 
redshift evolution of the rms mass fluctuation o-g inferred from Lya absorbers [20[ and from the power 
spectrum of density fluctuations measured from galaxies' peculiar velocities [21[. Current uncertainties 
are still too large to allow these measurements to discriminate among alternative cosmological scenarios, 
(e.g. [22, 23[). 

On-going redshift surveys like VIPERS [24[ or BOSS [25[ will certainly provide more stringent constraint 
and will be able to test those models that deviate most from the standard cosmological model. However, 
only next generation large-scale redshift surveys at z ss 1 and beyond like EUCLID [26[ or BigBOSS [27[ 
will provide an efficient way to discriminate competing dark energy models. 

The growth rate s clearly depends on the cosmological model. It has been found in several works 
[28-32[ that a simple yet effective parameterization of s captures the behavior of a large class of models. 
Putting 

s^^l, (3) 

where f2m(z) is the matter density in units of the critical density as a function of redshift, a value 
7 « 0.545 reproduces well the ACDM behavior while departures from this value characterize different 
models. For instance the DGP is well approximated by 7 « 0.68 [33, 34[ while viable models of f{R) 
are approximated by 7 « 0.4 for small scales and small redshifts [5, 6[. This simple parameterization 
is however not flexible enough to accommodate all cases. A constant 7 cannot for instance reproduce a 
growth rate larger than s = 1 in the past (as we have in f{R) and scalar-tensor models) allowing at the 
same time s < 1 at the present epoch if flm < 1. Even in standard cases, a better approximation requires 
a slowly- varying, but not strictly constant, 7. 

In addition, the measures of the growth factor obtained from redshift distortions require an estimate of 
the galaxy bias, which can be obtained either independently, using higher order statistics (e.g. [19, 35[) 
or inversion techniques [36[ , or self consistently, by assuming some reasonable form for the bias function 
a priori (for instance, that the bias is independent of scale, as we will assume here). 

The goal of this paper is to forecast the constraints that future observations can put on the growth 
rate. In particular we use representative assumptions for the parameters of the EUCLID survey to 
provide a baseline for future experiments and we focus on the following issues, i) We assess how well 
one can constrain the bias function from the analysis of the power spectrum itself and evaluate the 
impact that treating bias as a free parameter has on the estimates of the growth factor. We compare 
the results with those obtained under the more popular approach of fixing the bias factor (and its error) 
to some independently-determined value, ii) We estimate how errors depend on the parameterization of 
the growth factor and on the number and type of degrees of freedom in the analysis. Hi) We explicitly 
explore the case of coupling between dark energy and dark matter and assess the ability of measuring 
the coupling constant. 

We do this in the context of the Fisher Matrix analysis. This is a common approach that has been 
adopted in several recent works, some of which exploring the case of a EUCLID-like survey as we do. 
We want to stress here that this work is, in fact, complementary to those of the other authors. Unlike 
most of these works, here we do not try to optimize the parameter of the EUCLID survey in order to 
improve the constraints on the relevant parameters, as in [37[. Instead, we adopt a representative sets 
of parameters that describe the survey and derive the expected errors on the interesting quantities. In 
addition, unlike [38[ and [39[, we do not explicitly aim to study the correlation between the parameters 
that describe the geometry of the system and the growth parameters, although in our approach we also 
take into account the degeneracy between geometry and growth. Finally, the main results of this paper 
are largely complementary to the work of [40[ that perform a more systematic error analysis that does 
not cover the main issues of our work. 

Although, as we mentioned, in general s might depend on scale, we limit this paper to an exploration 
of time-dependent functions only. Forecasts for specific forms of scale-dependent growth factor motivated 
by scalar-tensor models are in progress and will be presented elsewhere. 

The layout of the paper is as follows. In the next section we will introduce the different parameteriza- 
tions adopted for the growth rate and for the equation of state of dark energy, together with the models 
assumed for the biasing function, and describe the different cosmological models we aim to discriminate. 
In sec. Ill we will briefly review the Fisher matrix method for the power spectrum and define the adopted 
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fiducial model. In sec. IV we will describe the characteristics of the galaxy surveys considered in this 
work, while in sec. V we will report our results on the forecast errors on the parameters of interest. 
Finally, in sec. VI we will draw our conclusions and discuss the results. 



II. MODELS 



The main scope of this work is to quantify the ability of future redshift surveys to constrain the 
growth rate of density fluctuations. In particular we want to quantify how this ability depends on the 
parameterization assumed for s and for the equation of state of the dark energy w and on the biasing 
parameter. For this reason we explore different scenarios detailed below. 



A. Equation of state 



• w -parameterization. In order to represent the evolution of the equation of state parameter w, we 
use the popular CPL parameterization [41, 42] 

w{z) = Wo + Wly— . (4) 

As a special case we will also consider the case of a constant w. 



B. Growth Rate 



As anticipated, in this work we assume that the growth rate, s, is a function of time but not of scale. 
Here we explore three different parameterizations of s: 

• s -parameterization. This is in fact a non-parametric model in which the growth rate itself is modeled 
as a step- wise function 5(2:) — Si, specified in different redshift bins. The errors are derived on Si 
in each i-th redshift bin of the survey. 

• ^-parameterization. As a second case we assume 

s EE n,^{zy^'^ . (5) 

where the j{z) function is parametrized as 

z 

7(z) = 70 +71— — . (6) 
z 

As shown by [43, 44[, this parameterization is more accurate than that of eq. (3) for both ACDM 
and DGP models. Furthermore, this parameterization is especially effective to distinguish between 
a wCDM model (i.e. a dark energy model with a constant equation of state) that has a negative 
71 (-0.020 < 71 < -0.016, for a present matter density 0.20 < ^1^,0 < 0.35) and a DGP model 
that instead, has a positive 71 (0.035 < 71 < 0.042). In addition, modified gravity models show a 
strongly evolving 7(2) [5, 43, 45[, in contrast with conventional Dark Energy models. As a special 
case we also consider 7 = constant (only when w also is assumed constant), to compare our results 
with those of previous works. 

• rj-parameterization. To explore models in which perturbations grow faster than in the ACDM case, 
like in the case of a coupling between dark energy and dark matter [4[, we consider a model in 
which 7 is constant and the growth rate varies as 

S^n^r.izni+Tl) , (7) 

where rj quantifies the strength of the coupling. The example of the coupled quintessence model 
worked out by [4[ illustrates this point. In that model, the numerical solution for the growth rate 
can be fitted by the formula (7), with rj = cj3'^, where /3c is the dark energy-dark matter coupling 
constant and best fit values 7 = 0.56 and c = 2.1. In this simple case, observational constraints 
over 77 can be readily transformed into constraints over P^. 
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C. Galaxy Biasing 

In the analysis of the redshift distortions, s{z) is degenerate with the bias function b{z). In absence 
of a well-established theory of galaxy formation and evolution, most analysis assume some arbitrary 
functional form for b{z). However, biasing needs to be neither deterministic nor linear. Stochasticity in 
galaxy biasing is supposed to have little impact on two-point statistics, at least on scales significantly 
larger than those involved with galaxy evolution processes [46]. On the other hand, deviations from 
linearity (which imply scale dependency) might not be negligible. Current observational constraints 
based on self consistent biasing estimators [19, 36[ show that nonlinear effects are of the order of a few 
to ~ 10 %, depending on the scale and the galaxy type [35, 47[. To account for current uncertainties in 
both modeling and measuring galaxy bias we consider the following choices for the functional form of b: 

• Redshift dependent bias. We assume b{z) = \/l + z (already used in [48[) since this function provides 
a good fit to line galaxies with luminosity Lh„ = 10*^ erg~^ s""'^ h~^ modeled by [49[ using 
the semi-analytic GALFORM models of [50[. We consider Hq, line objects since they will likely 
constitute the bulk of galaxies in the next generation slitless spectroscopic surveys like EUCLID. 
This Hjj luminosity roughly corresponding, at z = 1.5, to a limiting flux of /hq > 4 x 10^^^ 
erg cm~^s^^. 

• Constant bias. For the sake of comparison, we will also consider the case of constant b ~ 1 
corresponding to the rather unphysical case of a redshift-independent population of unbiased mass 
tracers. 



D. Reference Cosmological Models 



As it will be better explained in the next section, to perform the Fisher Matrix analysis we need to 
adopt a fiducial cosmological model. We choose the one recommended by the Dark Energy Task Force 
(DETF) [51[. In this "pseudo" ACDM model the growth rate values are obtained from eq. (3) with 
7 — 0.545 and f2„j(z) is given by the standard evolution 

a„(z) = a„,o(i + 2)^777^ ' (8) 

H (zj^ 

where (the subscript will generally denotes the present value) 



H{zr = Hi 



f7„,o(l + zY + VLk{l + zf + (1 - a„,o - f^fc) exp p / 1 + Wo + w 



l + z 1 



dz 



(9) 

Then ^m{z) is completely specified by setting i7m,o = 0.25, VLk — wq — —0.95, wi ~ 0. We wish 
to stress that regardless of the parameterization adopted, our fiducial cosmology is always chosen as 
the DETF one. In particular we choose as fiducial values 71 = and 77 = 0, when the corresponding 
parameterizations are employed. 

One of the goals of this work is to assess whether the analysis of the power spectrum in redshift- 
space can distinguish the fiducial model from alternative cosmologies, characterized by their own set of 
parameters (apart from Vlmfl which is set equal to 0.25 for all of them). The alternative models that we 
consider in this work are: 

• DGP model. We consider the flat space case studied in [52[. When we adopt this model then we 
set 7o — 0.663, 71 = 0.041 [43[ or 7 = 0.68 [33[ and w — —0.8 when 7 and w are assumed constant. 

• f{R) model. Here we consider the one proposed in [53[, depending on two parameters, n and A, 
which we fix to n = 2 and A = 3. In this case we assume 70 = 0.43, 71 = —0.2, values that apply 
quite generally in the limit of small scales (provided they are still linear, see [5[) or 7 = 0.4 and 
w = -0.99. 

• coupled dark energy (CDE) model. This is the coupled model proposed by [54, 55[. In this case we 
assume 70 = 0.56, rj = 0.056 (this value comes from putting /3c — 0.16 as coupling, which is of the 
order of the maximal value allowed by CMB constraints) [56[. As already explained, this model 
cannot be reproduced by a constant 7. 
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III. FISHER MATRIX ANALYSIS 



In order to constrain the parameters, we use the Fisher matrix method [57] (see [58] for a review), that 
we apply to the power spectrum analysis in redshift space following ]9]. For this purpose we need an 
analytic model of the power spectrum in redshift space as a function of the parameters that we wish to 
constrain. The analytic model is obtained in three steps, (i ) First of all we compute with CMBFAST ]59] 
the linear power spectrum of the matter in real space at z = 0, Por{k), choosing a reference cosmology 
where the parameters to be given as input (i.e. ilm.oh^, ^b.oh^, h, Ug also employed in the Fisher matrix 
analysis, plus the other standard parameters required by the CMBFAST code) are set to the values given 
in the III column of Tab. I while for the normalization of the spectrum we use erg = 0.8. (ii) Second, we 
model the linear redshift-space distortions as 

Pobs(2; k, m) = Por{k) ^f}'^^^^\ G'{z)b'{z){l + /3{z)f,Y + P,{Z) . (10) 

D^{z)Hf{z) 

Here the subscript F indicates quantities evaluated on the fiducial model. In this expression H{z) is 
the expansion history in Eq. 9, D{z) is the angular diameter distance, G{z) the growth factor and Ps{z) 
represents a scale-independent offset due to imperfect removal of shot-noise. Finally /3(z) is the redshift 
distortion parameter and the term (1 + f3fj?)'^ is the factor invoked by [60] to account for linear distortion 
in the distant-observer's approximation, where fi is the direction cosine of the wavenumber k with respect 
to the line of sight. As shown in ]10, 11] and recently in ]37], the inclusion of growth rate information 
reduces substantially the errors on the parameters, improving the figure of merits. (Hi ) As a third and 
final step we account for nonlinear effects. On scales larger than (^ 100/i^^Mpc) where we focus our 
analysis, nonlinear effects can be represented as a displacement field in Lagrangian space modeled by an 
elliptical Gaussian function. Therefore, following [61, 62], to model nonlinear effect we multiply Por(fc) 
by the factor 



exp ■ 



(11) 



where 'S± and represent the displacement across and along the line of sight, respectively. They are 
related to the growth factor G and to the growth rate s through = SqG and = SoG(l -I- s). 
The value of Sq is proportional to as- For our reference cosmology where erg = 0.8 [63], we have 
Eo = ll/i^^Mpc. 

The observed power spectrum in a given redshift bin depends therefore on a number of parameters, 
denoted collectively as pi, such as the Hubble constant at present h, the reduced matter and baryon 
fractions at present, f2m,o^^ and ftt.oh^, the curvature density parameter H.^, the spectral tilt Ug plus the 
parameters that enter in the parameterizations described in the previous section: wq, wi (or simply w); 
7o, 7i (or 7) and rj. They are listed in Tab. I and are referred to as "Cosmological parameters". These 
parameters will be left free to vary while we always fix as=0.8 since the overall amplitude is degenerate 
with growth rate and bias. The other free parameters depend on the redshift. They are listed in the 
lower part of Tab. I and include the expansion history H{z), the growth factor G{z), the angular diameter 
distance D{z), the shot noise Ps{z), the growth rate s{z), the redshift distortion parameter f3{z) and the 
galaxy bias b{z). 

Given the model power spectrum we calculate, numerically or analytically, the derivatives 



dlnPpbs 
dpi 



(12) 

F 



evaluated at the parameter values of the fiducial model and we obtain for the i-th redshift bin all the 
elements of the Fisher matrix through [64] 



Purvey 1 (-^^) 



''-%^^J_, V dp. dp, 

where 



V'eff(A:, p) 



nP{k,n) "'^ 
nP{k,p) + 1 
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is the effective volume of the survey sampled at the scale k along the direction fi. Vgurvcy and n represent 
the volume of the survey and the mean number density of galaxies in each redshift bin. 

As a fiducial model we assume a "pseudo" ACDM with wq = —0.95; the differences with the standard 
Wq — —1.0 ACDM model are rather small. For example, in the case of the j -parameterization, our fiducial 
model has 70 — 0.545, 71 = whereas the standard ACDM model has 70 = 0.556, 71 — 0.018 [43]. To 
summarize, our fiducial model is the same model recommended by the Dark Energy Task Force [51], i.e.: 
n^,^o = 0.25, ^ 0.0445, f7f 0, /i^ = 0.7, nf = 1, = -0.95, wf = 0, 7-^ = 0.545, Pf = 0. In 
addition, we assume that 7^ = 0, 77^ = 0. The fiducial values for the redshift dependent parameters are 
computed in every bin through the standard Friedmann-Robertson- Walker relations 



H^{z) 

mohP 


= n^^^,{i + zf + {i~nf^,,){i + zf('^<)\'" , 


(15) 


D^{z) 


X 1 r dz' 


(16) 


s^{z) 




(17) 


G^{z) 




(18) 


/3^(z) 


bP{z) ' 


(19) 


b^{z) 


= 1 or b^{z) = \/l + z , 


(20) 




= . 


(21) 



At this point our analysis is performed in two ways, according to the choice of z-dependent parameters 
that characterize the power spectrum: 

• Internal bias method. 

We assume some fiducial form for b{z) (z-dependent or constant) and express the growth function 
G{z) and the redshift distortion parameter /3(z) in terms of the growth rate s (see eqs. (22), (2)). 
When we compute the derivatives of the spectrum (eq. (12)), b{z) and s{z) are considered as 
independent parameters in each redshift bin. In this way we can compute the errors on b (and s) 
self consistently by marginalizing over all other parameters. 

• External bias method. 

In this case we also assume the same forms for b{z) as in the Internal bias case but we do not 
explicit G{z) and /3(z) in terms of s. The independent parameters are now the product G{z) ■ b(z) 
(if we considered them separately, the Fisher matrix would result singular) and /3(z). In this case 
we compute the errors over /3(z) marginalizing over all other parameters. Since we also marginalize 
over (G • 6)^, in this case we cannot estimate the error over b from the Fisher matrix. Thus, in 
order to obtain the error over s (related to /3 through s = /3 ■ b) with standard error propagation, 
we need to assume an "external" error for b{z). We allow the relative error Ab/b to be either 1% or 
10%, two values that bracket the ranges of expected errors contributed by model uncertainties and 
deviations from linear biasing. 



IV. MODELING THE REDSHIFT SURVEY 

The main goals of next generation redshift surveys will be to constrain the Dark Energy parameters 
and to explore models alternative to standard Einstein Gravity. For these purposes they will need 
to consider very large volumes that encompass z ~ 1, i.e. the epoch at which dark energy started 
dominating the energy budget, spanning a range of epochs large enough to provide a sufficient leverage 
to discriminate among competing models at different redshifts. The additional requirement is to observe 
some homogeneous class of objects that are common enough to allow a dense sampling of the underlying 
mass density field. 
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Cosmological parameters in Pobs{z; k, n) 

Reduced total matter density 

Reduced baryon density 

Curvature density 

Hubble constant at present 

Primordial fluctuation slope 

Constant growth index 

'y -parameterization parameters 

r] -parameterization parameters 




7 

70, 7i 
7, V 



Fiducial values 



0.25- (0.7)^ 
0.0445 ■ (0.7)2 



0.545, 
0.545, 





0.7 
1 

0.545 



Redshift dependent parameters 
Hubble parameter 
Angular diameter distance 
Growth rate 
Growth factor 

Redshift distortion parameter 

Shot noise 

Bias 



logH 
logD 
logs 
logG 
log/3 

Ps 

log 6 



1, ^TT^ 



eq. (15) 
eq. (16) 
eq. (17) 
eq. (18) 
eq. (19) 







Table I: Parameters. 



As anticipated in the introduction, in this paper we consider as a reference case the spectroscopic 
survey proposed by the EUCLID collaboration [26]. We stress that our aim is not to focus on this 
particular redshift survey and assess how the constraints on the relevant parameters depends on the survey 
characteristics in order to optimize future observational strategies. On the contrary, under the hypothesis 
that next-generation space-based all-sky redshift surveys will be similar to the EUCLID spectroscopic 
survey, we consider the latter as a reference case and estimate how the expected errors on the bias, growth 
rate, coupling constant and other relevant quantities will change when one consider slightly different 
observational setups. For this purpose we take advantage of the huge effort made by the EUCLID team 
to simulate the characteristic of the target objects and compute the expected selection function and 
detection efficiency of the survey and adopt the same survey parameters presented in [65] . 

Here we consider a survey covering a large fraction of the extragalactic sky (|&| > 20°), corresponding 
to - 20000 deg2 capable to measure a large number of galaxy redshifts out to z ~ 2. A promising 
observational strategy is to target emitters at near-infrared wavelengths (which implies z > 0.5) since 
they guarantee both relatively dense sampling (the space density of this population is expected to increase 
out to z ^ 2) and an efhcient method to measure the redshift of the object. The limiting flux of the 
survey should be the tradeoff between the requirement of minimizing the shot noise, the contamination 
by other lines (chiefly among them the [Oii[ line), and that of maximizing the so-called efficiency e, i.e. 
the fraction of successfully measured redshifts. To minimize shot noise one should obviously strive for a 
low ffux. Indeed, the authors in [65[ found that a limiting ffux /hq > 1 x 10~^^ erg cm~^s~"'^ would be 
able to balance shot noise and cosmic variance out to z — 1.5. However, simulated observations of mock 
Ha galaxy spectra have shown that e ranges between 30 % and 60% (depending on the redshift) for a 
limiting ffux /h„ > 3 x 10""'^^ erg cm~^s^^ [26[. Moreover, contamination from [Oii[ line drops from 12% 
to 1% when the limiting ffux increases from 1 x 10^^® to 5 x 10^^^ [65[. Taking all this into account we 
adopt a conservative choice and consider three different surveys characterized by a limiting ffux of 3, 4 
and 5 x 10~^^ ergcm~^s~^. 

We use the number density of H^ galaxies at a given redshift, n{z), estimated in [65[ using the latest 
empirical data and obtained by integrating the Hq luminosity function above the minimum luminosity 
set by the limiting ffux LHo,min. = 47rI?2^(z)^/H„ where -Dl(z) is the luminosity distance. To obtain 
the effective number density one has to account for the success rate in measuring galaxy redshifts from 
Ha lines. The effective number density is then obtained by multiplying n{z) by the already mentioned 
efficiency, e. In the range of redshifts and ffuxes considered in this work the value of e varies in the 
interval [30%, 50%[ (see Fig. A. 1.4 of [26[). 

In an attempt to bracket current uncertainties in modeling galaxy surveys, we consider the following 
choices for the survey parameters: 

• Reference case (ref.). Limiting ffux: /hq > 4 x 10~^® erg cm^^s^^, which gives the galaxy number 
densities listed in Col. 2 of Tab. II. The efficiency is set to e — 0.5. 
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z 


ni(2)xl0~^ 


n2(2)xl0"^ 


n3(2)xlO~^ 


0.5-0.7 


4.69 


3.56 


2.8 


0.7-0.9 


3.33 


2.42 


1.84 


0.9-1.1 


2.57 


1.81 


1.33 


1.1-1.3 


2.1 


1.44 


1.03 


1.3-1.5 


1.52 


0.99 


0.68 


1.5-1.7 


0.92 


0.55 


0.35 


1.7-1.9 


0.54 


0.29 


0.17 


1.9-2.1 


0.31 


0.15 


0.08 



Table II; Expected galaxy number densities in units of (/i/Mpc)^ for EUCLID survey. 

• Optimistic case (opt.). Limiting flux: /ho > 3 x 10^^^ erg cni^^s^^, which gives the galaxy number 
densities Hsted in CoL 1 of Tab. IL The efficiency is set to e = 0.5. 

• Pessimistic case (pess.). Limiting flux: /hq > 5 x 10~^^ erg cm~^s~"'^, which gives the galaxy 
number densities listed in Col. 3 of Tab. ff. The efficiency is set to e = 0.3. 

The total number of observed galaxies ranges from 3 • 10'^ (pess.) to 9 • lO"^ (opt.). For all cases we assume 
that the relative error on the measured redshift is az = 0.001, independent of the limiting flux of the 
survey. 

V. RESULTS 

In this section we present the main result of the Fisher matrix analysis that we split into two sections 
to stress the different emphasis given in the two approaches. We note that in all tables below we always 
quote errors at 68% probability level and draw in the plots the probability regions at 68% and/or 95% 
(denoted for shortness as 1 and 2a values). Moreover, in all figures, all the parameters that are not shown 
have been marginalized over or fixed to a fiducial value when so indicated. 

A. s-parameterization 

This analysis has two main goals: that of figuring out our ability to estimate the biasing parameter and 
that of estimating the growth rate with no assumptions on its redshift dependence. The total number 
of parameters that enter in the Fisher matrix analysis is 45: 5 parameters that describe the background 
cosmology (f2m,o^^, ^^fc,o^^, h, n, ^2fe) plus 5 z-dependent parameters specified in 8 redshift bins evenly 
spaced in the range z = [0.5,2.1]. They are Ps{z), D{z), H{z), s{z), b{z) in the internal bias case, while 
we have /3{z) and G{z) ■ b{z) in the place of s{z) and b{z) when we use the external bias method. 

The subsequent analysis depends on the bias method adopted. 

• In case of the internal bias method, the fiducial growth function G{z) in the (i + l)-th redshift bin 
is evaluated from a step- wise, constant growth rate s(z) as 

To obtain the errors on Si and bt we compute the elements of the Fisher matrix and marginalize 
over all other parameters. In this case one is able to obtain, self-consistently, the error on the bias 
and on the growth factor at different redshifts, as detailed in Tab. Ill and Tab. IV respectively. 

Tab. Ill illustrates one important result: through the analysis of the redshift-space galaxy power 
spectrum in a next-generation EUCLID-like survey, it will be possible to measure galaxy biasing 
in Az = 0.2 redshift bins with less than 3.5% error, provided that the bias function is independent 
on scale. This fact can be appreciated in Fig. 1 in which we show the expected relative error as 
a function of redshift for both b{z) functions and for the survey Pessimistic case. Errors are very 
similar in all but the outermost redshift shells. We show the Pessimistic case since with a more 
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0.015 



0.010 
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Z 

Figure 1: Relative errors over the bias parameter as a function of redshift, computed through the Fisher matrix 
analysis with the "internal bias" method for the pessimistic case. The blue solid line refers to the fiducial bias b = 1, 
while the purple dashed line refers to 6 = \/l + z. The precision in measuring the bias has a little dependence 
on the b{z) form: errors are very similar (the discrepancy is less than 1%) but in the outermost redshift shells 
(where however is less than 3%). 

favorable survey configuration, like the Reference case, the errors would be almost identical. In 
addition we find that the precision in measuring the bias has a little dependence on the b{z) form. 
The largest discrepancy between the b{z) = 1 and b{z) = ^/l + z cases is ^ 3% and refers to the 
expected errors on the growth rate at z = 2 in the Pessimistic case. Differences are typically much 
smaller for all other parameters or, for s{z) at lower redshifts or with a more favorable survey 
setup. Given the robustness of the results on the choice of b{z) in the following we only consider 
the b{z) = + z case. 

In Fig. 2 we show the errors on the growth rate s as a function of redshift, overplotted to our 
fiducial ACDM (green solid curve). The three sets of error bars are plotted in correspondence of 
the 8 redshift bins and refer (from left to right) to the Optimistic, Reference and Pessimistic cases, 
respectively. The other curves show the expected growth rate in three alternative cosmological 
models: fiat DGP (red dashed curve), f{R) (blue dotted curve) and CDE (purple, dot-dashed 
curve). This plot clearly illustrates the ability of next generation surveys to distinguish between 
alternative models, even in the less favorable choice of survey parameters. 

• In case of the external bias method we marginalize over the overall amplitude (G • 5)^. Since, 
in this case, we cannot find errors self-consistently, we assume that bias has been determined a 
priori with errors per redshift bin of 1% and 10%, two values that should bracket the expected 
range of uncertainties. We note that the external bias method can be considered more conservative, 
especially in the case of large errors although we see no obvious reason why it should be preferred to 
the internal bias method that seems to provide similar results. Indeed, the errors on s relative to the 
1% bias error listed in Table V are quite similar to those of the internal bias case. As expected, errors 
on s increase significantly when the bias is known with 10% accuracy rather than 1%. However, 
even in this case, one keeps the ability of distinguishing between most of the competing cosmological 
models at la level, as shown in Fig. 3. 

The main results of this section can be summarized as follows. 

1. The ability of measuring the biasing function is not too sensitive to the characteristic of the survey 
(b{z) can be constrained to within 1.5% in the Optimistic scenario and up to 3.5% in the Pessimistic 
one) provided that the bias function is independent on scale. Moreover, the precision in measuring 
the bias has a very little dependence on the b(z) form. 

2. The growth rate s can be estimated to within 1-3% in each bin for the Reference case survey with 
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b(z) 


— 1 (internal) 






Vl + 2 (internal) 


CTb 


z 


Oh 


ref. 


opt. 


pess. 




ref. 


opt. 


pess. 


0.014 


0.012 


0.017 


0.6 


0.012 


0.011 


0.015 


0.014 


0.012 


0.017 


0.8 


0.012 


0.011 


0.016 


0.013 


0.012 


0.017 


1.0 


0.012 


0.011 


0.015 


0.012 


0.011 


0.016 


1.2 


0.011 


0.010 


0.014 


0.012 


0.011 


0.015 


1.4 


0.011 


0.010 


0.013 


0.012 


0.010 


0.025 


1.6 


0.010 


0.010 


0.013 


0.012 


0.010 


0.019 


1.8 


0.010 


0.010 


0.014 


0.016 


0.011 


0.033 


2.0 


0.011 


0.010 


0.017 



Table III: la marginalized errors for the bias in each redshift bin obtained with the "internal bias" method. 



h{z) 


= 1 (internal) 






h{z) = 


Vl + ^ (internal) 
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ref. 


opt. 


pess. 






ref. 


opt. 


pess. 


0.011 


0.010 


0.013 
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0.73 


0.010 


0.010 


0.012 


0.011 


0.010 


0.015 


0.8 


0.78 
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0.013 


0.011 


0.016 


1.0 


0.83 


0.012 


0.010 


0.015 


0.013 


0.012 


0.018 


1.2 


0.86 


0.013 


0.011 


0.017 


0.014 


0.012 


0.019 


1.4 


0.89 


0.014 


0.012 


0.019 


0.016 


0.013 


0.023 


1.6 


0.91 


0.016 


0.013 


0.022 


0.019 


0.015 


0.032 


1.8 


0.92 


0.018 


0.015 


0.028 


0.027 


0.018 


0.057 


2.0 


0.93 


0.024 


0.017 


0.044 



Table IV: la marginalized errors for the growth rates in each redshift bin (Fig. 2) obtained with the "internal 
bias" method. 



no need of estimating the bias function b{z) from some dedicated, independent analysis using higher 
order statistics [19] or full-PDF analysis [36]. 

3. If the the bias were measured to within 1% in each slice, then the error over s would be very similar 
(just 1-2% larger) to that obtained by the internal estimate of h{z). 

4. The estimated errors on s depend weakly on the fiducial model of b{z). 



b{z) = 1 (external) 
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Ab/b = 10% 
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ref. 


opt. 
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ref. 
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ref. 


opt. 


pess. 


ref. 


opt. 


pess. 


0.014 


0.013 


0.015 


0.074 


0.074 


0.074 


0.6 


0.015 


0.014 


0.016 


0.074 


0.074 


0.074 
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0.012 


0.014 


0.079 


0.079 


0.079 


0.8 


0.013 


0.013 


0.015 


0.079 


0.079 


0.079 


0.012 


0.011 


0.014 


0.083 


0.083 


0.083 


1.0 


0.013 


0.012 


0.014 


0.083 


0.083 


0.083 


0.012 


0.012 


0.015 


0.086 


0.086 


0.087 


1.2 


0.013 


0.012 


0.015 


0.086 


0.087 


0.087 


0.014 


0.012 


0.018 


0.089 


0.089 


0.089 


1.4 


0.014 


0.013 


0.027 


0.089 


0.089 


0.089 


0.017 


0.014 


0.025 


0.092 


0.091 


0.094 


1.6 


0.016 


0.014 


0.022 


0.092 


0.091 


0.093 


0.023 


0.017 


0.042 


0.095 


0.093 


0.101 


1.8 


0.019 


0.016 


0.031 


0.094 


0.093 


0.097 


0.036 


0.023 


0.082 


0.099 


0.096 


0.124 


2.0 


0.027 


0.019 


0.052 


0.097 


0.095 


0.107 



Table V: lo marginalized errors for the growth rates in each redshift bin (Fig. 3) obtained with the "external 
bias" method. 
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Figure 2: Expected constraints on the growth rates in each redshift bin (using the "internal bias" method). The 
upper panel refers to fe = 1, while the lower panel to 6 = \/l + z. For each z the central error bars refer to 
the Reference case while those referring to the Optimistic and Pessimistic case have been shifted by -0.015 and 
+0.015 respectively. The growth rates for four different models are also plotted: ACDM (green solid curve), flat 
DGP (red dashed curve), f{R) model (blue dotted curve) and a model with coupling between dark energy and 
dark matter (purple, dot-dashed curve). In this case it will be possible to distinguish these models with next 
generation data. 



B. Other parameterizations. 



In this section we assess the ability of estimating s{z) when it is expressed in one of the parametrized 
forms described in Section II B. More specifically, we focus on the ability of determining 70 and 71, in the 
context of the ^ -parameterization and 7, 77 in the rj -parameterization. In both cases the Fisher matrix 
elements have been estimated by expressing the growth factor as 



G{z) = 5q exp 



(1-f 77) / 



dz' 
l + z' 



(23) 



where for the ^-parameterization we fix 77 = 0. In this section we adopt the internal bias approach and 
assume that h{z) = \/\ -\- z since, as we have checked, in the case of h{z) = 1 one obtains very similar 
results. 
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Figure 3: Expected constraints on the growth rates in each redshift bin (using the "external bias" method), 
assuming for the bias a relative error of 1% (upper panel) and 10% (lower panel). For each z the central error 
bars refer to the Reference case while those referring to the Optimistic and Pessimistic case have been shifted by 
-0.015 and +0.015 respectively. The growth rates for four different models are also plotted: ACDM (green solid 
curve), flat DGP (red dashed curve), f{R) model (blue dotted curve) and a model with coupling between dark 
energy and dark matter (purple, dot-dashed curve). Even in the case of large errors (10%) for the bias it will be 
possible to distinguish among three of these models with next generation data. 



• ^-parameterization. We start by considering the case of constant 7 and w in which we set 7 = 7^ = 
0.545 and w — — —0.95. As we will discuss in the next Section, this simple case will allow us to 
cross-check our results with those in the literature. In Fig. 4 we show the marginalized probability 
regions, at 1 and 2a levels, for 7 and w. The regions with different shades of green illustrates 
the Reference case for the survey whereas the blue long-dashed and the black short-dashed ellipses 
refer to the Optimistic and Pessimistic cases, respectively. Errors on 7 and w are listed in Tab. VI 
together with the corresponding figures of merit [FOM] defined to be the squared inverse of the 
Fisher matrix determinant and therefore equal to the inverse of the product of the errors in the 
pivot point, see [51]. Contours are centered on the fiducial model. The blue triangle and the blue 
square represent the flat DGP and the f{R) models' predictions, respectively. It is clear that, 
in the case of constant 7 and w, the measurement of the growth rate in a EUCLID-like survey 
will allow us to discriminate among these models. These results have been obtained by fixing the 
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case 






FOM 




ref. 


0.02 


0.02 


2115 




opt. 


0.019 


0.019 


2806 


Qk fixed 


pess. 


0.03 


0.03 


1296 



Table VI: Numerical values for la constraints on parameters in Fig. 1 and figures of merit. Here we have fixed 
Qk to its fiducial value. 

curvature to its fiducial value ilk — 0. If instead, we consider curvature as a free parameter and 
marginalize over, the errors on 7 and w increase significantly, as shown in Table VII, and yet the 
precision is good enough to distinguish the different models. For completeness, we also computed 
the fully marginalized errors over the other Cosmological parameters for the reference survey, given 
in Tab. VIII. 

As a second step we considered the case in which 7 and w evolve with redshift according to eqs. (6) 
and (4) and then we marginalize over the parameters 71, wi and il^.. The marginalized probability 
contours are shown in Fig. 5 in which we have shown the three survey setups in three different 
panels to avoid overcrowding. Dashed contours refer to the z-dependent parameterizations while 
red, continuous contours refer to the case of constant 7 and w obtained after marginalizing over 
fife. Allowing for time dependency increases the size of the confidence ellipses since the Fisher 
matrix analysis now accounts for the additional uncertainties in the extra-parameters 71 and wi; 
marginalized error values are in columns cr^^^^^ ^, f-mmarg 1 ^f Tab. IX. We note, however, that errors 
are still small enough to distinguish the fiducial model from the f{R) and DGP scenarios. 

We have also projected the marginalized ellipses for the parameters 70 and 71 and calculated 
their marginalized errors and figures of merit, which are reported in Tab. X. The corresponding 
uncertainties contours are shown in Fig. 6. Once again we overplot the expected values in the f{R) 
and DGP scenarios to stress the fact that one is expected to be able to distinguish among competing 
models, irrespective on the survey's precise characteristics. 

As a further test we have estimated how the errors on 70 depend on the number of parameters 
explicitly involved in the Fisher matrix analysis. Fig. 7 shows the expected la errors on 7 (Y- 
axis) as a function of the number of parameters that are fixed when computing the element of the 
Fisher matrix (the different combinations of the parameters are shown on the top of the histogram 
elements). We see that error estimates can decrease up to ~ 50 % when parameters are fixed to 
some fiducial value, or are determined independently. 

• rj -parameterization. 

We have repeated the same analysis as for the ^-parameterization taking into account the possibility 
of coupling between DE and DM i.e. we have modeled the growth factor according to eq. (7) and 
the dark energy equation of state as in eq. (4) and marginalized over all parameters, including ilj,. 
The marginalized errors are shown in columns cr-y^^^^^^i '^Wmarg 2 Tab. IX and the significance 
contours are shown in the three panels of Fig. 8 which is analogous to Fig. 5. The uncertainty 
ellipses are now larger than in the case of the "/-parameterization and show that DGP and f{R) 
models could be rejected at > la level only if the redshift survey parameter will be more favorable 
than in the Pessimistic case. 

Marginalizing over all other parameters we can compute the uncertainties in the 7 and rj parameters, 
as listed in Tab. XI. The relative confidence ellipses are shown in the left panel of Fig. 9. This plot 
shows that next generation EUCLID-like surveys will be able to distinguish the reference model 
with no coupling (central, red dot) to the CDE model proposed by [56] (white square) only at the 
1-1.5 a level. 

Finally, in order to explore the dependence on the number of parameters and to compare our results 
to previous works, we also draw the confidence ellipses for wq, wi with three different methods: i) fixing 
7o,7i to their fiducial values and marginalizing over all the other parameters; ii) marginalizing over all 
parameters plus 70,71 but fixing flk', Hi) marginalizing over all parameters but wq, wi. As one can see 
in Fig. 10 and Tab. XII this progressive increase in the number of marginalized parameters reflects in a 
widening of the ellipses with a consequent decrease in the figures of merit. These results are in agreement 
with those of other authors (e.g. [37, 40]). 
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Figure 4: 7-parameterization. 1 and 2a marginalized probability regions for constant 7 and w: the green (shaded) 
regions are relative to the Reference case, the blue long-dashed ellipses to the Optimistic case, while the black 
short-dashed ellipses are the probability regions for the Pessimistic case. The red dot marks the fiducial model; 
two alternative models are also indicated for comparison. 
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Figure 5: 7-parameterization. 1 and 2a marginalized probability regions obtained assuming constant 7 and w 
(red solid curves) or assuming the parameterizations (6) and (4) and marginalizing over 71 and wi (black dashed 
curves); marginalized error values are in columns a^^^^^ j, o"m,„apg 1 Tab. IX. Yellow dots represent the fiducial 
model, the triangles a f{R) model and the squares mark the flat DGP. 



The results obtained this Section can be summarized as follows. 

1. If both 7 and w are assumed to be constant and setting fifc = then, a redshift survey described 
by our Reference case will be able to constrain these parameters to within 4% and 2%, respectively. 

2. Marginalizing over ftk degrades these constraints to 5.5% and 4% respectively. 

3. If w and 7 are considered redshift-dependent and parametrized according to eqs (6) and (4) then 
the errors on 70 and Wq obtained after marginalizing over 71 and Wi increase by a factor 4, 5, 
i.e. we expect to measure 70 and wq with a precision of 13-15% and 11-14% respectively, where 
the interval reflects the uncertainties in the characteristic of the survey. With this precision we 
will be able to distinguish the fiducial model from the DGP and f{R) scenarios with more than 2a 
significance. 
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bias 


case 


0"-, 




FOM 




ref. 


0.03 


0.04 


1179 


b = Vl + z 


opt. 


0.03 


0.03 


1568 




pess. 


0.04 


0.05 


706 



Table VII: Numerical values for la constraints on parameters 7 and w (no parameterization), relative to the red 
ellipses in Figs 5, 8 and figures of merit. Here we have marginalized over fi^. 





case 












b = y/l + z 


ref. 


0.024 


0.008 


0.002 


0.01 


0.02 



Table VIII: Numerical values for marginalized la constraints on Cosmological parameters using constant 7 and 
w. 

4. The ability to discriminate these models with a significance above 2(t is confirmed by the confidence 
contours drawn in the 70-71 plane, obtained after marginalizing over all other parameters. 

5. If we allow for a coupling between dark matter and dark energy, and we marginalize over rj rather 
than over 71, then the errors on 70 and wq are almost identical to those obtained in the case of the 
"f -parameterization. However, our ability in separating the fiducial model from the CDE model is 
significantly hampered: the confidence contours plotted in the 7-77 plane show that discrimination 
can only be performed wit 1-1. 5cr significance. 



VI. CONCLUSIONS 

In this paper we addressed the problem of determining the growth rate of density fiuctuations from 
the estimate of the galaxy power spectrum at different epochs in future redshift survey. As a reference 
case we have considered the proposed EUCLID spectroscopic survey modeled according to the latest, 
publicly available survey characteristics [26, 65]. In this work we focused on a few issues that we regard 
as very relevant and that were not treated in previous, analogous Fisher Matrix analysis mainly aimed at 
optimizing the survey setup and the observational strategy. These issues are: i) the ability in measuring 
self-consistently galaxy bias with no external information and the impact of treating the bias as an 
extra free parameter on the error budget; ii) the impact of choosing a particular parameterization in 
determining the growth rate and in distinguishing dark energy models with very different physical origins 
(in particular we focus on the ACDM, f{R) and the DGP, models that are still degenerate with respect 
to present growth rate data) ; Hi) the estimate of how errors on the growth rate depend on the degrees 
of freedom in the Fisher matrix analysis; iv) the ability of estimating a possible coupling between dark 
matter and dark energy. 

The main results of the analysis were already listed in the previous Section, here we recall the most 
relevant ones. 

1. With the "internal bias" method we were able to estimate bias with 1% accuracy in a self consistent 
way using only galaxy positions in redshift-space. The precision in measuring the bias has a very 
little dependence on the functional form assumed for b(z). Measuring b with 1% accuracy will be 
a remarkable result also from an astrophysical point of view, since it will provide a strong, indirect 
constraint on the models of galaxy evolution. 



bias 


case 






FOM 




°"'"'marg.2 


FOM 




ref. 


0.08 


0.11 


241.1 


0.09 


0.11 


104.1 




opt. 


0.07 


0.09 


323.5 


0.07 


0.09 


138.2 




pess. 


0.11 


0.14 


142.5 


0.11 


0.15 


61.6 



Table IX: la marginalized errors for parameters 7 and w expressed through 7 and rj parameterizations. Columns 
■yo,margi,wo,margi refer to marginalizatiou over 71, Wi (Fig. 5) while columns jo,7narg2,wo,marg2 refer to marginal- 
ization over 1^,101 (Fig. 8). 
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71 




Figure 6: 7-parameterization. 1 and 2a marginalized probability regions for the parameters 70 and 71, relative 
to the Reference case (shaded yellow regions), to the Optimistic case (green long-dashed ellipses), and to the 
Pessimistic case (black dotted ellipses). Red dots represent the fiducial model, blue squares mark the DGP while 
triangles stand for the f{R) model. Then, in the case of 7-parameterization, one could distinguish these three 
models (at 95% probability). 



bias 


case 


(^70 




FOM 




ref. 


0.08 


0.17 


178.4 




opt. 


0.07 


0.15 


227.5 




pess. 


0.11 


0.22 


112.4 



Table X: Numerical values for la constraints on parameters in Fig. 6 and figures of merit. 



2. We have demonstrated that measuring the amplitude and the slope of the power spectrum in 
different z-bin allows to constrain the growth rate with good accuracy, with no need to assume an 
external error for h{z). In particular, we found that s can be constrained at Icr to within 3% in each 
of the 8 redshift bin from z = 0.5 to 2.1. This result is robust to the choice of the biasing function 
b{z). The accuracy in the measured s will be good enough to discriminate among the most popular 
competing models of dark energy and modified gravity. 

3. Taking into account the possibility of a coupling between dark matter and dark energy has the 
effect of loosening the constraints on the relevant parameters, decreasing the statistical significance 
in distinguishing models (from ^ 2cr to < l.Scr). Yet, this is still a remarkable improvement over the 
present situation, as can be appreciated from Fig. 9 where we compare the constraints expected by 
next generation data to the present ones. Moreover, the Reference survey will be able to constrain 
the parameter ry to within 0.04. Reminding that we can write 77 = 2.1/3^ [4], this means that the 
coupling parameter /3c between dark energy and dark matter can be constrained to within 0.14, 
solely employing the growth rate information. This is comparable to existing constraints from the 
CMB but is complementary since obviously it is obtained at much smaller redshifts. A variable 



bias 


case 


a-, 


an 


FOM 




ref. 


0.08 


0.04 


464.1 


b = ^/l + z 


opt. 


0.07 


0.03 


608.2 




pess. 


0.11 


0.05 


280.3 



Table XI: Numerical values for la constraints on parameters in Fig. 9 and figures of merit. 
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Figure 7: The bars represent the errors on the growth index 70 obtained using the "/-parameterization and fixing 
an increasing number of cosmological parameters as indicated over each bar and marginalizing over the others. 
The progressive increase in the number of fixed parameters reflects in a decrease of the error. 
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Figure 8: r;-parameterization. 1 and 2a marginalized probability regions obtained assuming constant 7 and w 
(red soHd curves) or assuming the parameterizations (7) and (4) and marginahzing over r] and wi (black dashed 
curves); marginalized error values are in columns a-y^^^^ 21 ""^marg 2 of Tab. X. Yellow dots represent the flducial 
model, the triangles stand for a f{R) model and the squares mark the flat DGP. 



coupling could therefore be detected by comparing the redshift survey results with the CMB ones. 

It is worth pointing out that, whenever we have performed statistical tests similar to those already dis- 
cussed by other authors in the context of a EUCLID-like survey, we did find consistent results. Examples 
of this are the values of FOM and errors for wq, wi, similar to those in [37, 40] and the errors on constant 7 
and w [40]. However, let us notice that all these values strictly depend on the parameterizations adopted 
and on the numbers of parameters fixed or marginalized over. In particular, we also found that all these 
constraints can be improved if one uses additional information from e.g. CMB and other observations. 









FOM 


70,71 fixed 


0.04 


0.24 


166 


fixed and marginalization over both 70,71 


0.12 


0.36 


97.3 


marginalization over all other parameters 


0.12 


0.43 


41.3 



Table XII: 1 a marginalized errors for the parameters Wq and Wi, obtained with three different methods (reference 
case, see Fig. 10 ). 



18 



T] 0.0 





-0.3 



7 



Figure 9: 7)-parameterization. Left panel: 1 and 2a marginalized probability regions for the parameters 7 and rj 
in eq. (7) relative to the reference case (shaded blue regions), to the optimistic case (yellow long-dashed ellipses) 
and to the pessimistic case (black short-dashed ellipses). The red dot marks the fiducial model while the square 
represents the coupling model. Right panel: present constraints on 7 and rj computed through a full likelihood 
method (here the red dot marks the likelihood peak) [4]; long-dashed contours are obtained assuming a prior for 



We made a first step in this direction in Fig. (7), which shows how the errors on a constant 7 decrease 
when progressively more parameters are fixed by external priors. 
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Figure 10: Errors on the equation of state. 1 and 2cr marginalized probability regions for the parameters wq and 
wi, relative to the reference case and constant bias 6=1. The blue dashed ellipses are obtained fixing 70,71 to 
their fiducial values and marginalizing over all the other parameters; for the red shaded ellipses instead, we also 
marginalize over 70,71 but we fix £7^ = 0. Finally, the black dotted ellipses are obtained marginalizing over all 
parameters but wq and wi. The progressive increase in the number of parameters reflects in a widening of the 
ellipses with a consequent decrease in the figures of merit (see Tab. XII). 
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